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Using numerical simulations we examine the static and dynamic properties of the recently proposed 
vortex liquid crystal state. We confirm the existence of a smectic-A phase in the absence of pinning. 
Quenched disorder can induce a smectic state even at T = 0. When an external drive is applied, a 
variety of anisotropic dynamical flow states with distinct voltage signatures occur, including elastic 
depinning in the hard direction and plastic depinning in the easy direction. We disuses the impli- 
cations of the anisotropic transport for other systems which exhibit depinning phenomena, such as 
stripes and electron liquid crystals. 
PACS numbers: 74.25.Qt 



Recently, a new state of vortex matter termed a vor- 
tex liquid crystal was proposed to occur in superconduc- 
tors with anisotropic vortex- vortex interactions [1]. In 
such systems, the vortex lattice first melts in the soft di- 
rection, giving rise to an intermediate vortex smectic-A 
state, followed at higher temperatures by a melting into 
a nematic state. The initial theoretical calculations of 
this transition combined an elastic model with the Lin- 
demann criterion for melting; however, the validity of this 
approach has been called into question [2,3]. The Lindc- 
mann criterion does not take into account the prolifera- 
tion of dislocations that is likely to occur in a smectic-A 
state, so a numerical investigation would be very useful 
both to determine whether the smectic-A state can occur 
in the vortex system as well as to examine the dynam- 
ics of vortex liquid crystals in the presence of disorder. 
The physics of the vortex liquid crystal state should be 
generic to the class of problems which can effectively be 
modeled as a two-dimensional (2D) system of particles 
with anisotropic repulsive interactions. Such a system 
has already been physically realized for magnetic col- 
loidal particles, where a smectic-A state was observed 
along with dislocations that have preferentially aligned 
Burgers vectors [4]. There is also considerable interest 
in electron liquid crystal states, which may arise when 
anisotropic interactions in classical electron crystals give 
rise to smectic and nematic states [5,6]. Evidence for 
such states has been observed in transport measurements 
which show hard and soft directions for flow [7,8]. 

Previous studies of vortex smectic states considered 
vortices interacting with some form of an underlying ID 
periodically modulated substrate [9] . A very similar sys- 
tem in which smectic states have been observed is col- 
loidal particles interacting with ID periodic substrates 
[10]. Each of these systems melts into an intermediate 
smectic- C state. When the underlying substrate is dis- 
ordered rather than periodic, application of an external 
drive induces an anisotropic fluctuating force that or- 
ganizes the vortices into a moving smectic state, where 
the dislocations in the vortex lattice are aligned with the 
direction of the applied drive [11,12]. In the case of the 
proposed vortex liquid crystal state, the anisotropy arises 



when the vortex cross section becomes elliptical due to an 
anisotropic superfluid stiffness which leads to different ef- 
fective masses in the three crystalline directions [1] . The 
theoretical calculations in Ref. [1] were performed for a 
system with no quenched disorder; however, real super- 
conductors often contain significant amounts of random 
pinning. It would be desirable to understand the trans- 
port properties of vortex liquid crystals in order to iden- 
tify signatures of the liquid crystal phase and seek new 
types of dynamical phenomena in these systems. Under- 
standing how quenched disorder affects an anisotropic 
system of repulsively interacting particles is also relevant 
to dynamics in electron liquid crystal states. 

To address these issues, we consider a 2D system of N v 
interacting vortices with periodic boundary conditions 
in the x and y directions. The overdamped equation of 
motion for a single vortex i is 

^ = fr + ff + ff + ff a) 

The damping constant 77 is set to unity. The vortex- 
vortex interaction force is f™ = Y^fli A v Ki{r.ij / \)r\j, 
where K\ is the modified Bessel function, which decays 
exponentially for large distances, A is the London pene- 
tration depth, A v is the vortex interaction prefactor, and 
Tij is the distance between vortices i and j. The Bessel 
function is appropriate for stiff, 3D vortex lines. We have 
also considered 1 jr interaction potentials appropriate for 
thin film superconductors as well as Yukawa interaction 
potentials for colloidal particles and find the same qual- 
itative features. The thermal force f[ arises from ran- 
dom Langevin kicks with the properties < ff >= and 
< f[{t)ff{t?) >= 2rjk B T5{t - t')Sij. The quenched dis- 
order ff is modeled as random pinning sites in the form 
of attractive parabolic traps of radius r p = 0.2 A and 
strength f p . The Lorentz driving force from an exter- 
nal applied current is f d . The system size is measured in 
units of A and the forces in terms of A v . The anisotropic 
interactions are introduced by multiplying the vortex- 
vortex interaction force in the x and y directions by a 
vector (C x ,Cy), where the anisotropy C = C x /C y . In 
this work we concentrate on the case C = 1/ \/l0 which 
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FIG. 1. (a,c,e) Black dots: vortices; black lines: vortex 
trajectories. (b,d,f) Delaunay triangulation, with topological 
defects (5 and 7-fold coordinated particles) marked as filled 
circles. (a,b) T = 0.5; (c,d) T = 1.2; (e,f) T = 1.35. 

is the value considered in Ref. [1]. We take the x axis to 
be the soft direction and the y axis as the hard direction. 

We first consider the case where the pinning and the 
external driving force are absent. In Fig. 1 we illustrate 
the melting of a 24A x 24A system with a vortex density 
of p v — 1.2/A 2 . Figure 1(a) shows the vortex positions 
(dots) and trajectories (lines) for a fixed period of time 
with a fixed T = 0.5, and Fig. 1(b) shows a corresponding 
Delaunay triangulation. At this temperature, the system 
remains in a crystalline state with no dislocations. The 
vortices are undergoing larger random displacements in 
the soft (x) direction than in the hard (y) direction; how- 
ever, there is no long time diffusion of the particles. Fig- 
ures 1(c) and 1(d) present the smectic-A state at T = 1.2. 
Here the trajectories have a ID liquid structure with 
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FIG. 2. The average particle displacements in each direc- 
tion, d x and d y , normalized by the lattice constant a, vs time, 
measured in molecular dynamics steps, (a) The smectic-A 
state at T = 1.21. (b) The nematic phase at T = 1.35. 

motion along the soft x direction and no significant trans- 
lation of the vortices in the y direction. The Delau- 
nay triangulation indicates the presence of dislocations 
with aligned Burgers vectors, which is characteristic of 
the smectic-A state. Figures 1(e) and 1(f) illustrate the 
vortex liquid phase at T = 1.35. The vortex trajecto- 
ries show clear diffusion in both the x and y directions, 
with more pronounced motion in the x direction. The 
dislocations are no longer aligned in a single direction, 
indicating the loss of long-range order in both the x and 
y directions. These results confirm that a smectic-A state 
can occur in a system of vortices with anisotropic inter- 
actions, as predicted by theoretical calculations [1]. We 
note that when the anisotropy ratio C is too small, the 
two-step melting transition illustrated here is lost. 

To further characterize the smectic state, in Fig. 2 we 
plot the average particle displacements for the x and 
y directions, d x =< X^f" \ x i(t) ~ x i(t')\ > /N v and 
d y =< X^f" \Vi{t)~ Ui{t')\ > /N v . In the smectic phase at 
T = 1.21, shown in Fig. 2(a), d x /a increases much more 
rapidly than d y /a, and does not saturate but increases 
to a value over 1, indicating that the vortices can diffuse 
more than a lattice constant in the x direction over time. 
This is due to the formation of dislocations which allow 
adjacent rows of vortices to slip past each other while 
remaining confined in the y direction. We note that the 
saturation value ofd y /a is approximately 1/5, larger than 
the Lindcmann criterion value of 1/10. Excess motion in 
the y direction occurs during a sliding event when two 
rows slip past each other and the vortices in each row 
are temporarily displaced in the direction perpendicular 
to the slip plane. This transverse motion is not large 
enough to permit the formation of dislocations aligned 
in the hard direction. In the nematic phase, shown in 
Fig. 2(b) at T = 1.35, d x still increases more rapidly 
than d y ; however, the continuous increase of both quan- 
tities indicates that the particles are diffusing throughout 
the entire system. 

We next consider the effect of random disorder by 
adding N p = 2N V randomly located pinning sites to the 
same system studied in Fig. 1, and then conducting a 
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FIG. 3. Regions in which the smectic and disordered phases 
occur in a system with quenched disorder as a function of 
temperature T and pinning strength f p . Dashed line roughly 
indicates the weak pinning region in which a 2D anisotropic 
Bragg glass forms. Inset: the density of six-fold coordinated 
particles Pa vs T for (top curve) f v = 0.025 and (bottom 
curve) f p = 0.2. 

series of simulations at varied T and varied pinning 
strength f p . For high temperatures we always obtain 
a nematic phase, while at T = for large f p we find a 
pinned nematic phase. At low T and low f p we observe a 
phase very similar to that shown in Fig. l(c,d), where the 
vortex lattice is oriented in the soft direction and there 
are a small number of aligned dislocations. We term this 
a pinned smectic-A phase. In Fig. 3 we indicate the re- 
gions in which the smectic and disordered phases appear 
as a function of temperature and pinning strength. The 
phase boundary is identified via the density of sixfold co- 
ordinated particles, Ps; the defect density is given by 1- 
Pq. In the crystal phase, there are no defects and P& = 1. 
In the smectic phase, P§ — 0.91 to 0.95, and in the ne- 
matic phase Pq > 0.8. In the inset of Fig. 3 we plot Pq 
vs T for two different disorder strengths. For f p = 0.025 
(upper line) the system is in the pinned smectic state at 
T = 0. As T increases, there is a clear transition to the 
disordered state, as indicated by the drop in in P 6 near 
T = 1.19. The lower line shows Pq for f p = 0.2, when the 
pinning is strong enough to disorder the system even at 
T = 0. These results suggest that weak random disorder 
can increase the extent of the regions where the smectic- 
A phase occurs when there are anisotropic interactions, 
by suppressing the crystalline phase at low temperatures 
and raising the melting temperature of the smectic state. 

It has been shown that dislocations are always present 
for weak disorder in two dimensional isotropic systems; 
however, the distance between the dislocations can be 
arbitrarily large compared to the range of the 




FIG. 4. (a) Pq vs driving force /<j for a system with 
f p — 0.04 at T = 0. Upper curve: Pg for id = fay- Lower 
curve: Pg for fd = /die. (b) Average velocities vs fd for the 
same system. Upper curve: V y for f d = fay- Lower curve: V x 
for f d — /dX. (c) V x and V y vs f d for a system with f p = 0.2. 
(d) A blowup of (c) in the region near depinning showing the 
crossing of the velocity force curves. 

translational order, so that for a wide range of tempera- 
tures and disorder strengths the system behaves as a 2D 
Bragg glass [14]. In the vortex liquid crystal case there 
are two length scales associated with the hard and soft di- 
rections, and the disorder induced dislocations first form 
in the soft direction. It may be possible that, on very 
large length scales, dislocations in the hard direction will 
also appear. For the parameters considered here, disloca- 
tions are present except at the lowest pinning strengths, 
roughly indicated by a dashed line in Fig. 3, where a 2D 
anisotropic Bragg glass forms. 

We next consider dynamical effects in the presence of 
pinning. In the smectic-A state, there should be dis- 
tinct transport signatures for the hard and soft direc- 
tions. When the disorder is strong enough to destroy the 
smectic phase, there may still be an anisotropic transport 
signature if the system retains some form of nematic or- 
der. We first consider the pinned smectic-A state found 
at f p = 0.04 and T = 0. We perform separate simulations 
for driving in the soft direction, = /<pc, and the hard 
direction, id = fd$, increasing the applied drive very 
slowly to avoid any transient effects. In Fig. 4(b) we plot 
V y = (1/N V ) < v v y > (upper curve) for driving in 

the hard direction and V x = (1/N V ) < Vx > (l° wer 
curve) for driving in the soft direction. Here, V x < V yi 
indicating that motion in the soft direction is easier ex- 
cept at very high drives when the effects of the pinning 
are washed out and the two curves come together. In 
Fig. 4(a) we plot Pq for the two different driving direc- 
tions. At fd = 0, Pe is slightly less than one due to the 
presence of a small number of dislocations in the smectic 
state. For depinning in the hard direction, P^ (upper 
curve), the system depins elastically without a prolifcra- 
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tion of defects, and the vortices do not exchange neigh- 
bors as they move. For driving in the easy direction, Pq 
(lower curve) drops substantially when the vortices depin 
plastically, and a portion of the vortices remain pinned 
while others flow past. The effective pinning is known to 
be higher for a soft system where defects can proliferate 
than for an elastic system [15]. A similar proliferation of 
defects has been associated with the so called peak effect, 
where there is a sudden increase in the effective pinning 
force as a function of temperature or applied magnetic 
field [15]. At high drive, the system shifts from plastic 
flow to a dynamically reordered state [12] as indicated by 
the increase in Pg in Fig. 4(a), as well as by the merging 
of V x and V y in Fig. 4(b). These results imply that, in 
the pinned smectic-A state, the depinning is elastic in the 
hard direction and plastic in the soft direction. Further, 
the elastic and plastic depinning transitions produce dif- 
ferent scaling responses in the velocity force curves. At 
depinning, the velocity scales with the driving force in 
the form V = (fa — fcY [13]. For the plastic flow regime 
we find (3 > 1.0 while in the elastic flow regime we find 
/3 < 1.0, in agreement with theoretical expectations. 

At finite temperatures and for f p large enough that 
we observe only plastic depinning in both directions, we 
observe that the critical depinning force in the soft direc- 
tion, is lower than the depinning force in the hard di- 
rection, even though V x < V y at intermediate drives. 
This implies that the anisotropic flow exhibits a reversal 
from V x > V y to V x < V y at low drives. We explicitly 
demonstrate this effect for a system with f p — 0.25 and 
T = 0.25 in Fig. 4(c,d). Here, the depinning is plastic in 
both directions, and the depinning forces are f£ = 0.015 
and f}! = 0.11. There are fewer dislocations for = fdf 
and the system reorders at f% = 0.2. For = /^x, the 
system does not reorder until f% = 0.8. There is a clear 
crossing of the velocity force curves at fd = 0.13 so that 
the flow is easier in the soft direction for fd < 0.13 and 
easier in the hard direction for fd > 0.13. In Fig. 4(d) we 
show a blowup of this region. The crossing of the velocity 
force curves can be understood by considering that the 
depinning in the soft direction is plastic. At low drives, 
individual vortices can be thermally activated, giving rise 
to creep. For driving in the hard direction, the depinning 
is elastic and individual vortex hopping is not possible, 
so that only collective creep can occur. In the case of 
the nematic phase, when there is some plastic flow in the 
y-direction, there is still a large correlated length scale 
that must move so thermal effects are greatly reduced. 
Thus, creep in the pinned smectic phase and pinned ne- 
matic phase is enhanced in the soft direction compared 
to the hard direction. Even for very high values of f p we 
observe an intermediate anisotropic response, suggesting 
the system can be considered a pinned nematic phase. 

In conclusion, we have performed simulations of the 
recently proposed vortex liquid crystal state where the 
vortex-vortex interactions are anisotropic. We find that, 



in the absence of disorder, the system shows an interme- 
diate melting into a smectic-A state as proposed theoret- 
ically in Ref. [1]. The smectic-A state contains a small 
fraction of dislocations which are all aligned in the soft 
direction. In the presence of disorder, a pinned smectic- 
A state can occur, and the system depins plastically in 
the soft direction but elastically in the hard direction. We 
predict that, for equal intermediate drives, the velocity in 
the soft direction will be lower than in the hard direction. 
At finite temperatures the creep is much higher in the soft 
direction due to the fact that individual vortex hopping 
can occur, whereas creep is suppressed in the hard direc- 
tion since the vortex motion is much more correlated. For 
high temperatures and disorder strengths, the system is 
disordered. For strong disorder the anisotropic transport 
should still be observable in the pinned nematic phase. 
We note that many of our results may apply to electron 
liquid crystals as well. 
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